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Abstract 

The influence of exponential magnetic field decay (MFD) on the spin evolution of 
isolated neutron stars is studied. The ROSAT observations of several X-ray sources, 
which can be accreting old isolated neutron stars, are used to constrain the expo- 
nential and power-law decay parameters. We show that for the exponential decay 
the ranges of minimum value of magnetic moment, //&, and the characteristic decay 
time, t d , ~ 10 29 - 5 > fj. b > 10 28 G cm 3 , ~ 10 8 > td > 10 7 yrs are excluded assum- 
ing the standard initial magnetic moment, (Mi = 10 30 G cm 3 . For these parameters, 
neutron stars would never reach the stage of accretion from the interstellar medium 
even for a low space velocity of the stars and a high density of the ambient plasma. 
The range of excluded parameters increases for lower values of (Iq. 

We also show, that, contrary to exponential MFD, no significant restrictions 
can be made for the parameters of power-law decay from the statistics of isolated 
neutron star candidates in ROSAT observations. 

Isolated neutron stars with constant magnetic fields and initial values of them 
less than (jlq ~ 10 29 G cm 3 never come to the stage of accretion. 

We briefly discuss the fate of old magnetars with and without MFD, and describe 
parameters of old accreting magnetars. 
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1 Introduction 



Astrophysical manifestations of neutron stars (NSs) are determined by 
their periods and magnetic fields. Four main evolutionary stages of iso- 
lated NSs can be singled out (see e.g. (Lipunov, 1992) for more details): 
the ejector, the propeller, the accretor and the georotator. 

On average NSs should have high spatial velocities due to an additional 
kick obtained during the supernova explosion (Lyne and Lorimer, 1994, 
Lorimer et al., 1997). The interstellar medium (ISM) accretion rate for high 
velocity objects should be rather low. However, recent population synthesis 
calculations (Popov et al., 2000) indicate that several old accreting NSs can 
be observed in the solar vicinity even for the space velocity distribution 
similar to one derived from radio pulsar observations. 

Magnetic field decay (MFD) in NSs is a matter of controversy. Many 
models of the MFD have been proposed starting from the first simple 
models (Gunn and Ostriker, 1970) up to the recent calculations (Sang and 
Chanmugam, 1990; Urpin and Muslimov, 1992). Observations of radio 
pulsars (Lyne et al., 1998) give no evidence for MFD with characteristic 
time scales td shorter than ~ 10 7 yrs. Here we suggest to use old accreting 
isolated neutron stars as probes of the models of field decay and try to put 
some limits on the parameters of the exponential and power-law MFD on 
a longer time scale assuming that some X-ray sources observed by ROSAT 
are indeed old accreting isolated NSs (Haberl et al., 1998; Neiihauser and 
Triimper, 1999; Schwope et al., 1999). 

This presentation is based mainly on two our recent papers (Popov and 
Prokhorov, 2000a,b). 

2 Calculations and results 

The main idea is to calculate the ejector time, i.e. a time interval spent by 
a NS on the ejector stage, for different parameters of the MFD and, using 
standard assumptions for the initial NS parameters, to compare this time 
with the Hubble time, tn- 

The ejector time, monotonically increases with increasing velocity 
of NS, v, and decreasing density of the ISM, n. For a constant magnetic 
field of a NS this relation takes the simple form: 
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tp(// = const) ~ 10 /i 30 n ' Vioyrs. (1) 

Using a high mean ISM density n ~ lcm~ 3 and a low space velocity of 
NSs (about the sound speed in the ISM), v ~ lOkms -1 , we arrive at the 
lower limit of tp. After the ejection stage has been over, the NS passes 
to the propeller stage and only after that can become an accreting X-ray 
source. The duration of the propeller stage tp is poorly known, but for a 
constant magnetic field is always less than tp, (see Lipunov and Popov, 
1995). Therefore if for some parameters of a NS tp exceeds the Hubble 
time tp ~ 10 10 yrs, it can not come to the accretion stage and hence can 
not underly the ROSAT INS candidate. 

We note, that if the initial magnetic moment of a NS is about fiQ ~ 10 29 
G cm 3 or smaller, than (in the case of constant field) this star never leave 
the ejector stage even for low velocity and high ISM density! So, significant 
part of INS for any velocity distribution can't become accretors at all. 

In addition, we assumed that NSs are born with sufficiently small rota- 
tional periods, po, and all have the same parameters of the MFD. We shall 
consider different initial surface magnetic field values. 

2.1 Exponential decay 

The field decay in this subsection is assumed to have an exponential shape: 

fi = fio ■ e~ t/td , for ji> jib (2) 

where is the initial magnetic moment (fi = ^B p R% s , here B p is the polar 
magnetic field, R^s is the NS radius), td is the characteristic time scale of 
the decay, and fib is the bottom value of the magnetic moment which is 
reached at the time t cr : 

= t d ■ In ( — ] , (3) 

and does not change after that. 

In Fig. 1 we show as an illustration the evolutionary tracks of NSs on 
P — B diagram for v = lOkms -1 and n = 1 cm~ 3 . Tracks start at t = 
when p = 20ms and fi = 10 30 Gcm 3 and end at t = tp = 10 10 yrs (for 
td = 10 7 yrs and td = 10 8 yrs) or at the moment when p = Pe (for 
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td = 10 9 yrs, td = 10 10 yrs and for a constant magnetic field) . The line 
with diamonds shows p = pe{B). 

The ejector stage ends when the critical ejector period, p E , is reached: 



n _ 11 r 1/2 -1/41/2 

p E — ll.o /i 30 n ^io s ? 



(4) 



where i>io = Jv* + *u 2 /10 kins" 1 . i> p is the NS's spatial velocity, u s and n are 
the sound velocity and density of the ISM, respectively. In the estimates 
below we shall assume v = 10 kms^ 1 and n = 1 cm" 3 . 

The initial NSs' spin periods should be taken much smaller than pe- 
Here to calculate duration of the ejection stage we assume po = s. To 
compute this time we used the magnetodipole formula: 



dp _ 2 47T 2 // 2 
dt 3 £>ic 3 
where \i can be a function of time. 

After a simple algebra we arrive at the following expression for 

-t d • In 
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where the coefficient T (which is just ts for constant magnetic field) is 
determined by the formula: 



T = 



3/c 



10 17 / 4 5/io>ro 1/2 Mn 1/2 s. 



(7) 



2fi V2vM 

Here M can be formally determined according to the Bondi equation for 
the mass accretion rate even if the NS is not at the accretion stage: 



M ~ 10 11 nv 10 3 gs~ 1 . 



(8) 



The results of calculations of for several values of fiQ and td are shown 
in Fig. 2. The right end points of all curves are limited by the values 
fib = fiQ. These points correspond to the evolution of a NS with constant 
magnetic field (see eq. (0)) and for them ts = T. One can see the increase 
of t E for evolution with a constant field for smaller initial fields. 
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If fib is small enough, the NS field has no time to reach the bottom value. 
In this case ts is determined by the 1st branch of equation (0) and does 
not depend on fib. In Fig. 2 this situation corresponds to the left horizontal 
parts of the curves. At 



the situation changes so that ts starts to depend on fib- In this region 
two counter-acting factors operates. On the one hand, the NS braking 
becomes slower with decreasing \x (see eq. (0)). On the other hand, the 
end period of the ejection pe becomes shorter Since ts < T at the left 
hand side horizontal part and (cZTe/g?/z&) | < 0, the right hand side of the 
curve must have a maximum. The first factor plays the main role to the 
right of the maximum. The magnetic field there rapidly falls down to fib 
at p <C Pe and most time NS evolves with the minimum field fi = /if, (this 
time period increases with decreasing fib). To the left of the maxium but 
before the horizontal part the NS's magnetic field reaches fi = fib with the 
spin period close to pe (the smaller fib, the closer) and soon after t = t cr 
the NS leaves the ejection stage. 

As it is seen from this Figure, for some combination of parameters ts is 
longer than the Hubble time. It means that such NSs never evolve further 
than the ejection stage. 

We argue that since accreting isolated NSs are really observed, the com- 
binations of td and fib for which no accreting isolated NS appear can be 
excluded for the progenitors of ROSAT X-ray sources. The regions of 
excluded parameters are plotted in Figs. 3 and 4. 

The hatched regions correspond to parameters for which ts is longer 
than 10 10 yrs, so a NS with such parameters never comes to the accretor 
stage and hence can not appear as an accreting X-ray source. In view of 
the fact of observations of accreting old isolated NSs by ROSAT satellite, 
this region can be called "forbidden" for a given ^o- 

In the "forbidden" region in Fig. 3, which is plotted for fiQ = 10 30 G cm 3 , 
all NSs reach the bottom field in a Hubble time or faster, and the evolution 
on late stages proceeds with the minimal field. The left hand side of the 
forbidden region is determined approximately by the condition 




PE(l^b) = p(t = t cr ). 



(9) 
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The right hand side of the region is roughly determined by the value 
of /ifr, with which a NS can reach the ejection stage for any i.e. this 
fib corresponds to the minimum value of with which a NS reaches the 
ejection stage without MFD. 

In Fig. 3 we also show the "forbidden" region for = 0.5 • 10 30 Gcm 3 
(dotted line). The dashed line in Fig. 3 shows that for all interesting 
parameters a NS with fiQ = 10 30 Gcm 3 reaches fib in less than 10 10 yrs. 
The dash-dotted line shows the same for fiQ = 0.5 • 10 30 Gcm 3 . The solid 
line corresponds to PE(f*>b) = p{t = t cr ), where t cr = td ■ ln(/io//i&). The 
physical sense of this line can be described in the following way. This line 
divides two regions: in the upper left region td are relatively long and fib 
relatively low, so NS can't reach bottom field during ejector stage; in the 
lower right region td are short and fib relatively high, so NS reach fib at the 
stage of ejection. 

Fig. 4 is plotted for fi = 10 29 Gcm 3 . For long t d (> 4 • 10 9 yrs) the NS 
cannot leave the ejection stage for any fib < fio- That's why in the upper 
part of the figure a horizontal "forbidden" region appears. 



2.2 Power- law decay 

Power-law (as also exponential) MFD is a widely discussed variant of NSs' 
field evolution. Power-law is a good fit for several different calculations of 
the field evolution (Goldreich & Reisenegger, 1992; Geppert et al., 2000). 
The power-law MFD can be described with the following simple formula: 

^ = -aB 1+a . (10) 
dt v 1 

So, we have two parameters of decay: a and a. As far as this decay is 

relatively slow for the most interesting values of a greater/about 1 (we 

use the same units as in (Colpi et al., 2000)), we don't specify any bottom 

magnetic field, contrary to what we made for more rapid exponential decay 

(Popov & Prokhorov, 2000a). Even for the Model C from (Colpi et al., 

2000) (see Table 1) with relatively fast MFD the magnetic field can decrease 

only down to ~ 10 8 G in 10 10 yrs (see Fig. 5). But for very small a 

the magnetic field can decay significantly during the Hubble time for any 

reasonable value of a. And, probably, it is useful to introduce in the later 

case a bottom field. 
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At the stage of ejection an INS is spinning down according to the mag- 
netodipole formula: PP « bB 2 . Here 6 = 3, values of magnetic field, B, 
Boq and Bq, are taken in units 10 13 G and time, t, in units 10 6 yrs (as in 
Colpi et al, 2000). 

In the table we show parameters of the Models A, B, C from (Colpi et 
al., 2000). Boo is the magnetic field calculated for t = t Hubble = 10 10 yrs and 
for the initial field Bq = 10 12 G. Models A and B correspond to ambipolar 
diffusion in the irrotational and the solenoidal modes respectively. Model 
C describes MFD in the case of the Hall cascade (models are valid mostly 
for relatively high values of magnetic field) . 



Table 1: Models A,B,C from (Colpi et al., 2000) 



Model 


A 


B 


C 


a 


0.01 


0.15 


10 


a 


5/4 


5/4 


1 


Boo 


« 1.9- 10 11 G 


« 2.4- 10 10 G 


« 10 8 G 



In Fig. 6 we show dependence of the ejector period, pe, and the asymp- 
totic period, poo, on the parameter a for a = 1 for different values of the 
initial magnetic field, Bq: 



Pe 



= 25.75V 2 n-V% 1 / 2 s , (11) 



Pi = ^-Bl a . (12) 
2 — aa 

Here pe was calculated for t = tnubbie = 10 10 yrs, i.e. for the moment, 
when B = Boo. 

It is clear from Fig. 6, that for the initial field greater/about 10 11 G low 
velocity INSs are able to come to the stage of accretion: for Bq = 10 11 G 
lines for poo and pe for the lowest possible velocity, 10 km/s, coincides. 

For power-law decay we can also plot "forbidden" regions on the plane 
a-a, where an INS for a given velocity for sure cannot come to the stage 
of accretion in the Hubble time (see Popov & Prokhorov, 2000b). If one 
also takes into account the stage of propeller (between ejector and accretor 
stages) it becomes clear, that "forbidden" regions for an INS which cannot 
reach the stage of accretion are even larger. 
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For the most interesting cases (Models A, B, C from (Colpi et al., 2000)) 
and v < 200 km/s INSs can reach the stage of accretion. It is an important 
point, that fraction of low velocity NSs is very small (Popov et al., 2000) 
and most of NSs have velocities about 200 km/s. 

3 Evolved magnetars 

In the last several years a new class of objects - highly magnetized NSs, 
"magnetars" (Duncan and Thompson, 1992) - became very popular in con- 
nection with soft 7-repeaters (SGR) and anomalous X-ray pulsars (AXP) 
(see Mereghetti and Stella, 1995; Kouveliotou et al., 1999; Mereghetti, 1999 
and recent theoretical works Alpar, 1999; Marsden et al., 2000; Perna et 
al., 2000). 

Magnetars come to the propeller stage with periods ~ 10 - 100 s in the 
Models A, B, C (see Fig. 2 in Colpi et al., 2000). Then their periods quickly 
increase, and NSs come to the stage of accretion with significantly longer 
periods, and at that stage they evolve to a so-called equilibrium period 
(Lipunov and Popov, 1995; Konenkov and Popov, 1997) due to accretion 
of the turbulent ISM: 

Peq ~ 2800Sr4 / V 2 /\4 3 ^C 2/3 M l :f 3 b (13) 

Here v t is a characteristic turbulent velocity, I - moment of inertia, M - 
INS's mass. This formula underestimate the period for relatively high v t , 
and relatively low v, because it assumes, that all external angular moment 
can be accreted by a INS. 

Isolated accretor can be observed both with positive and negative sign of 
p (Lipunov and Popov, 1995). Spin periods of INSs can differ significantly 
from p eq contrary to NSs in disc-fed binaries, and similar to NSs in wide 
binaries, where accreted matter is captured from giant's stellar wind. It 
happens because spin-up/spin-down moments are relatively small. 

As the field is decaying the equilibrium period is decreasing, coming to 
~ 28 sec when the field is equal to 10 10 G (we note here recently discovered 
objects RX J0420.0-5022 (Haberl et al, 2000) with spin period - 22.7 s). 

It is important to discuss the possibility, that evolved magnetar can 
appear also as a georotator. It happens if: 
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v > 300B u 1/5 n 1/w km/s. (14) 

For all values of a and a that we used NSs at the end of their evolution 
(t = 10 10 yrs) have magnetic fields < 10 12 G for wide range of initial fields, 
so they never appear as georotators if v < 480 km/s for n = 1cm -3 . But 
without MFD magnetars with B > 10 15 G and velocities v > 100 km/s 
can appear as georotators. 

Popov et al. (2000) showed, that georotator is a rare stage for INSs, 
because an INS can come to the georotator stage only from the propeller 
or accretor stage, but all these phases require relatively low velocities, and 
high velocity INSs spend most of their lives as ejectors. This situation is 
opposite to binary systems, where a lot of georotators are expected for fast 
stellar winds (wind velocity can be much faster than INS's velocity relative 
to ISM). 

Without MFD magnetars also can appear as accreting sources. In that 
case they can have very long periods and very narrow accretion columns 
(that means high temperature). Such sources are not observed now. Ab- 
sence of some specific sources associated with evolved magnetars (binary 
or isolated) can put some limits on their number and properties (dr. V. 
Gvaramadze drew our attention to this point). 

At the accretion part of INSs' evolution periods stay relatively close to 
p eq (but can fluctuate around this value), and INSs' magnetic fields decay 
down to ~ 10 10 — 10 11 G in several billion years for the Models A and 
B. It corresponds to the polar cap radius about 0.15 km and temperature 
about 250 - 260 eV (the same temperature, of course, can appear for INSs 
evolving with constant field), higher than for the observed INS candidates 
with temperature about 50 - 80 eV. We calculate the polar cap radius, 
Reap = Rns\J{Rns/Ra) {Ra - Alfven radius), with the following formula: 

R cap = 6 • 10 3 Sr 3 2/ V/ Vo 3/7 ^s fl cm - ( 15 ) 

The temperature can be even larger, than it follows from the formula above 
as far as for very high field matter can be channeled in a narrow ring, so 
the area of the emitting region will be just a fraction of the total polar cap 
area. 

As the field is decreasing the radius of the polar cap is increasing, and the 
temperature is falling. Sources with such properties (temperature about 
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250-260 eV) are not observed yet (Schwope et al., 1999). But if the num- 
ber of magnetars is significant (about 10% of all NSs) accreting evolved 
magnetars can be found in the near future, as far as now we know about 
5 accreting INS candidates (Treves et al., 2000; Neiihauser and Trumper, 
1999), and their number can be increased in future, p measurements are 
necessary to understand the nature of such sources, if they are observed. 

Recently discovered object RX J0420.0-5022 (Haberl et al, 2000) with 
the spin period ~ 22.7 s, can be an example of an INS with decayed 
magnetic field accreting from the ISM, as previously RX J0720. 4-3125. Due 
to relatively low temperature, 57 eV, its progenitor cannot be a magnetar 
for power-law MFD (Models A,B,C) or similar sets of parameters, because a 
very large polar cap is needed, which is difficult to obtain in these models. 
Of course RX J0420. 0-5022 can be explained also as a cooling NS. The 
question "are the observed candidates cooling or accreting objects?" is 
still open (see Treves et al., 2000). If one finds an object with p > 100 
s and temperature about 50 - 70 eV it can be a strong argument for its 
accretion nature, as far as such long periods for magnetars can be reached 
only for very high initial magnetic fields for reasonable models of MFD and 
other parameters. 

4 Discussion and conclusions 

We tried to evaluate the region of parameters which are excluded for models 
of the exponential and power-law MFD in NSs using the fact of observations 
of old accreting isolated NSs in X-rays. 

For the exponential decay the intermediate values of td (~ 10 7 — 10 8 yrs) 
in combination with the intermediate values of fit (~ 10 28 - 10 29 - 5 G cm ) 
for fiQ = 10 30 Gcm 3 can be excluded for progenitors of isolated accreting 
NSs because NSs with such parameters would always remain on the ejector 
stage and never pass to the accretion stage. 

For high NSs should reach ts even for td < 10 8 yrs. For weaker fields 
the "forbidden" region becomes wider. The results are dependent on the 
initial magnetic field ^o, the ISM density n, and NS velocity v. 

In fact the limits obtained are very strong, because we did not take into 
account that NSs can spend some significant time (in the case of MFD) 
at the propeller stage (the spin-down rate at this stage is very uncertain, 
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see the list of formulae, for example, in (Lipunov and Popov, 1995) or 
(Lipunov, 1992)). 

Note that there is another reason for which a very fast decay down to 
small values of fib can also be excluded, as far as this would lead to a huge 
amount of accreting isolated NSs in drastic contrast with observations. 
This situation is similar to the "turn-off" of the magnetic field of a NS 
(i.e., quenching any magnetospheric effect on the accreting matter). So for 
any velocity and density distributions we should expect significantly more 
accreting isolated NS than we know from ROSAT observations (of course, 
for high velocities X-ray sources will be very dim, but close NSs can be 
observed even for velocities ~ 100 km s _1 ). 

For power-law MFD (contrary to exponential decay) we cannot put 
serious limits on the parameters of decay with the ROSAT observations of 
INS candidates as far as for all plausible models of power-law MFD INSs 
from low velocity tail are able to become accretors. For more detailed 
conclusions a NS census for power-law MFD is necessary, similar to non- 
decaying and exponential cases (Popov et al., 2000). 

An interesting possibility of observing evolved accreting magnetars ap- 
pear both for the case of MFD and for constant field evolution. These 
sources should be different from typical present day INS candidates ob- 
served by ROSAT. Existence or absence of old accreting magnetars is very 
important for the whole NS astrophysics. 

We conclude that the existence of several old isolated accreting NSs 
observed by ROSAT (if it is the correct interpretation of observations), 
can put important bounds on the models of the MFD for isolated NSs 
for exponential decay (without influence of accretion, which can stimulate 
field decay). These models should explain the fact of observations of ~ 10 
accreting isolated NSs in the solar vicinity. Here we can not fully discuss 
the relations between decay parameters and X-ray observations of isolated 
NSs without detailed calculations. What we showed is that this connection 
should be taken into account and made some illustrations of it, and future 
investigations in that field would be desirable. 
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Figure 1: Tracks on P-B diagram. Tracks are plotted for bottom polar magnetic field 
8 • 10 10 G, initial polar field 2 • f0 12 G, NS velocity lOkms -1 , ISM density lcm" 3 and 
different tj. The last point of tracks with different corresponds to the following NS 
ages: 1 10 yrs for = 10 7 and = 10 8 yrs; 1.5 x 10 9 yrs for = 10 9 yrs; ~ 2 • 10 9 yrs for 
td = 10 10 yrs. The initial period is assumed to be po — 0.020 s. The line with diamonds 
shows the ejector period, pe- 
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Figure 2: Ejector time £g (in billion years) vs. the bottom value of the magnetic moment. 
The curves are shown for two values of the initial magnetic moment: 10 30 Gcm 3 (upper 
curves) and 10 31 Gcm 3 . 
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Figure 3: The characteristic time scale of the MFD, t^, vs. bottom magnetic moment, 
/lb. In the hatched region t# is greater than 10 10 yrs. The dashed line corresponds to 
tH — td • In (a«o / A*&) j where tjj = 10 10 years. The solid line corresponds to PE(Hb) = 
pit = t cr ), where t cr = td ■ ln(/z //ib). Both the lines and hatched region are plotted 
for /io = 10 30 Gcm~ 3 . The dash-dotted line is the same as the dashed one, but for 
fXo — 5 • 10 29 Gcm 3 . The dotted line shows the border of the "forbidden" region for 
fio = 5 • 10 29 Gcm 3 . 
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Figure 4: The characteristic time scale of the MFD, vs. bottom magnetic moment, 
yUb- In the hatched region is greater than 10 10 yrs. The dashed line corresponds to 
tu = t^-ln (/Iq/ fib), where tu = 10 10 yrs. The solid line corresponds to PE^b) = p(t = t cr ), 
where t cr = t d ■ In {fi Q / fi h ). Both lines and region are plotted for /x = 10 29 Gem -3 . 




Figure 5: Power-law MFD. Model A: a = 0.01, a = 1.25; solid line with circles. Model B: 
a = 0.15, a = 1.25; dashed line with squares. Model C: a = 10, a = 1; long-dashed line 
with diamonds. Models were described in details in Colpi et al. (2000) (see also Table 1). 
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Figure 6: Periods vs. parameter a for different values of the initial magnetic field: 
10 11 , 10 12 , 10 13 , 10 14 G. 
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